Decoherence of Bell states by local interactions with a dynamic spin environment 
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We study the evolution of a system of two qubits, each of which interacts locally with a spin chain 
with nontrivial internal Hamiltonian. We present a new exact solution to this problem and analyze 
the dependence of decoherence on the distance between the interaction sites. In the strong coupling 
regime we find that decoherence increases with increasing distance. In the weak coupling regime 
the dependence of decoherence with distance is not generic (i.e., it varies according to the initial 
state). Decoherence becomes independent of distance when the latter is over a saturation length I. 
Numerical results for the Ising chain suggest that the saturation scale is related to the correlation 
length ^. For strong coupling we display evidence of the existence of non-Markovian effects (such as 
environment-induced interactions between the qubits). As a consequence the system can undergo a 
quasiperiodic sequence of "sudden deaths and revivals" of entanglement, with a time scale related 
to the distance between qubits. 
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I. INTRODUCTION 

Understanding the process of decoherence [l|, 0, Q is 
not only important from a fundamental point of view but 
also essential to design good error correction strategies to 
prevent the collapse of quantum computers [3| . When a 
composite system interacts with an environment, deco- 
herence typically generates loss of entanglement between 
fragments. This process may depend on many details 
such as the nature of the system-environment interac- 
tion, the existence of non-trivial spatial correlations in 
the environment, the internal environmental dynamics, 
etc. In this paper we make a step towards understanding 
decoherence when different subsystems are locally cou- 
pled to a common environment with non-trivial internal 
dynamics and spatial correlations. 

In a previous publication Q we analyzed the decoher- 
ence induced on a single qubit by the interaction with an 
environment formed by a spin chain with an XY Hamil- 
tonian (see [51] for references on decoherence by spin en- 
vironments). Here, we consider a system of two qubits, 
each of which interacts with a different site of an XY spin 
chain. We solve this model generalizing previous results 
and obtain exact expressions for the evolution of the re- 
duced density matrix of the two qubits for a large family 
of initial states (that includes any of the four, maximally 
entangled. Bell states). Our work goes beyond previous 
studies extending the usual "central qubit" model, where 
each qubit is homogeneously coupled to all the sites of 
the chain In fact, our results enable us to analyze 
the way in which decoherence and disentanglement de- 
pend on the distance between the interaction sites. In 
this context we will study two limiting cases: when both 
qubits interact with the same site in the chain, and when 
the interaction sites are widely separated. More inter- 
estingly, we will investigate intermediate situations and 
examine the nature of the transition between the two 
limits. We will show that the environment correlation 
length plays a crucial role in setting the typical length 



for which decoherence stops depending on distance. We 
will also analyze the way in which the evolution of the 
system depends on the strength of the coupling between 
the system and the environment. As we will see, weak 
and strong coupling regimes are drastically different. For 
strong coupling we will show that the chain provides a 
medium through which the qubits can coherently inter- 
act. In such a case, the entanglement between them 
may exhibit quasiperiodic events of "sudden deaths" and 
"sudden revivals" 0] with a time scale that depends on 
the distance between qubits. 

The paper is organized as follows: in Section |TT] we in- 
troduce the model, defining the Hamiltonians for the sys- 
tem, the environment, and the coupling between them. 
We also present the main formulas we will use to deter- 
mine the decay of quantum coherence. In Sections [Hi] and 
IIVI we study the loss of entanglement between the qubits 
as a consequence of their interactions with the environ- 
ment, in the cases of weak and strong coupling with the 
chain respectively. Finally, in Section |V] we summarize 
our results. 



II. THE MODEL AND ITS ANALYTIC 
SOLUTION 

We study the decoherence induced on a system of two 
spin 1/2 particles (the qubits) by the coupling to an envi- 
ronment formed by a chain of spin 1/2 particles (Fig. 
H]). We neglect the self- Hamiltonian of the system. The 
total Hamiltonian is H = He + Hint, where He is the 
Hamiltonian of an XY spin chain: 
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(1) 

Here Xj , Yj , Zj denote the three Pauli operators acting 
on the j-th site of the chain, and we assume periodic 
boundary conditions. The parameter 7 determines the 
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FIG. 1: A system of two non-interacting qubits {A and B), 
suffering decoherence from the local interaction with distant 
sites of an environment chain. 



anisotropy in the x^y plane and A gives a magnetic field 
in z direction (7 = 1 corresponds to the Ising chain with 
transverse field). This model is critical for 7 with 
|A| < 1, and for A = ±1. 

Each of the two qubits in the system interacts locally 
with a certain spin of the chain: qubit A interacts with 
spin 0, and qubit B with spin d. The interaction Hamil- 
tonian is: 

H^nt = -5(|1)(1U^0 + (2) 

with |0) and |1) the two eigenstates of the Z Pauli op- 
erator. This has a simple interpretation: the coupling 
to the qubits induces a change in the effective magnetic 
field at sites and d. Thus, if the system is in state \ab) 
(a, b € {0, 1}) the environment evolves with an effective 
Hamiltonian Hab given by: 

Hab^Hc-g{aZo + bZd). (3) 

We also assume that the initial state of the "universe" 
formed by system and environment is of the form: 

PSEiO)=po®\Eo){Eo\ (4) 

where the initial state of the environment \Eq) is the 
ground state of the effective Hamiltonian Hqq [g]. 

Our goal is to study the evolution of p, the reduced 
density matrix of the two-qubit system (obtained from 
the state of the universe by tracing out the environment). 
Because of the special form of the Hamiltonian, the tem- 
poral dependence of p can be formally obtained as fol- 
lows. In the basis of eigenstates of Za and Zb, p can be 
written as: 

P{i)^ H PabM{1^)\ab){cd\. (5) 

abed— 0,1 

The evolution of the matrix elements of p is given by: 

Pab,cd{t) = Pab,cd{^){E^W"^''e-''''^''\Eo) (6) 

(we take h — 1). As states in the basis commute with 
the total Hamiltonian, the diagonal terms Pab,ab remain 
constant. On the other hand, each off-diagonal term in 
the reduced density matrix is modified by a factor of 
absolute value between and 1. This factor corresponds 



to the overlap between two different evolutions of the spin 
chain. Following we denote the square modulus 

of this factor as the Loschmidt echo: 

LabMt) = |(i?o|e'^^^*e-'^°^*|£;o)|', (7) 

which obviously satisfies Lat.cd — Lcd,ab- 

Below, we will show how to compute the echoes £00, ab 
and Loi,io- These two echoes are enough to obtain the 
reduction of the off-diagonal elements in the reduced den- 
sity matrix when the initial state is one of the four Bell 
states (or any state satisfying pu^i = pii.io = 0). The 
computation of each of these two echoes is slightly dif- 
ferent. In fact, to obtain £oo,afc we first note that one of 
the evolution operators in ([7]) acts trivially. Therefore, 
the echo is equal to the survival probability of the initial 
state after being evolved with the Hamiltonian Hab, i.e. 

Loo,ab{t) = \{Eo\e-'"-^'\Eo)\' (8) 

In turn, the echo Loi.io can be calculated noticing that 
the corresponding effective Hamiltonians are related by 
a translation: 

Hoi = T^'HioT-^, (9) 

where T is the one-site translation operator in the chain. 
As the ground state of the Hamiltonian Hqq is an eigen- 
statc of this operator, the echo can be written as: 

Loi.iQit) = |(i?o|e'^^''*T'^e"'^-*|i?o)|' (10) 

To obtain the echoes we will make use of the fact that 
the full quantum evolution for each Hamiltonian Hat can 
be exactly solved. This can be done by mapping the 
Hamiltonians Hab of the chain onto a fermion system by 
means of the Jordan- Wigner transformation |ll| : 




Z, - 2c]cj-l. (13) 



Using this, up to a correction term associated to bound- 
ary effects, the Hamiltonians can be written as: 

Hab = - XI + i + c]_^iCj) + 7(c]c] + i + Cj + iCj) + 

j 

+Af)(2c]c,-1)} (14) 

with A^"'''' = A -f g{a5jfi -f bSj^d) (the extension to qubits 
interacting with more than one site is trivial) . 

The Hamiltonians Hab depend quadratically on the an- 
nihilation and creation operators. Therefore they can 
be diagonalized by linear (Bogoliubov) transformations 
defining new creation and annihilation operators which 
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we will denote as 7y('''')^ jytCafe)^ Furthermore, as all these 
transformations are linear, the operators corresponding 
to different values of the labels (ab) can also be connected 
by Bogoliubov transformations. 

The echoes we want to compute can be written in terms 
of the matrices involved in these Bogoliubov transforma- 
tions. For example, in Appendix A it is shown that: 
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(t) = det(g + he'^') 



(15) 



Here A is a diagonal Nx N matrix containing the energies 
of the normal modes of the Hamiltonian Hat- g and h 
correspond to the transformation connecting the particles 
that diagonalize the unperturbed Hamiltonian Hqo and 
the effective Hamiltonian Hab- 
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ab) 



(16) 



This equation is useful since it expresses the echo as the 
determinant of an iV x matrix, which can be efficiently 
computed (the number of operations is polynomial in A'^). 
This formula is a new version of the one used in [9'] , where 
the Loschmidt echo was written in terms of the two-point 
correlators of the environment chain (i.e., as the deter- 
minant of a 2 A^ X 2 A^ matrix) . 

The way to compute the echo Loi.io is shown in Ap- 
pendix B, using the fact that the translation operator 
T is Gaussian in the fermion operators. The result is, 
once again, the determinant of an A^ x A^ matrix, though 
somewhat more complicated; 



Loi,io(t)= det(g'r-V^+h'T'*h'^^ 



(17) 



Here g', h' are time-dependent complex matrices related 
to Bogoliubov transformations between different sets of 
particle operators, and r is the diagonal matrix expres- 
sion of the translation operator T. 

Decoherence for initial states of the form \lp) — ajOO) 4- 
/3|11) or \ip) = a\Ql) -|-/3|10) is described by Loo,ii and 
Loi.io respectively. In any case, the relevant echo L(t) for 
the Bell-like state \if) will determine not only the process 
of purity decay but also the way in which the two qubits 
become disentangled. In fact, if one considers an initial 
mixture of the form po = pW) {v\ + (1 ~p) 1/4, the purity 
of the state is the following function of the echo L{t): 



+ p^[l - 2\a(3\\l - L)]. (18) 



(here, L is either Loi,io or Loo,ii depending on \ip)). The 
entanglement between the two qubits can be measured by 
the negativity Af obtained from the sum of the negative 
eigenvalues of the partial transpose of p (A/" 9^ is an 
indicator of non-separability of the two-qubit state) [l^] . 
For the mixed initial state proposed: 



Af{t) = max <^ 0, p\af3\ ^/ijj) - 



1 



P 



(19) 



For pure initial states {p = 1), the state is entangled 
whenever L ^ 0; in general, the state becomes disentan- 
gled when the echo is below a threshold which depends on 
p. In this way the phenomenon of entanglement sudden 
death (ESD) may take place (3]. 

We will analyze the dependence of the echo on the 
distance between the interaction sites. Two simple lim- 
its exist. First, if both qubits couple to the same site 
{d = 0) then the echo Lgi^io becomes trivially 1 as the 
two effective Hamiltonians are identical. On the other 
hand, Loo, 11 is simply the echo of a single qubit inter- 
acting with one site with twice the interaction strength 
(a case studied in [§]). The long distance limit is also 
easy to understand: if the chain is sufficiently large we 
expect its effect to be equal to the one obtained when 
each qubit interacts with an independent environment. 
Then, both Loo,ii and Loi,io approach Lq ^ (where Lo,i 
is the single-qubit echo) . This long distance regime exists 
only if the back-action of the qubits on the environment 
is small. For strong back-action, one expects the environ- 
ment to induce effective interactions between the qubits. 
As we will see, the long distance limit will be identifi- 
able easily in the weak coupling case, whereas the strong 
coupling regime will be plagued with effects associated 
to environment-induced interactions. 



III. RESULTS: WEAK COUPLING 

In this Section we study the case when the coupling g is 
small compared to the interaction strength between chain 
sites, showing results for the reference value g = 0.1. 
Fig. [2]shows the echo ioo,ii as a function of time for the 
Ising chain (7 = 1) and for several values of the trans- 
verse field A and the distance d. For A < 1, the echo has 
small amplitude oscillations (with frequency of order 1) 
about a value which is roughly independent of distance. 
For A > 1 decoherence is an order of magnitude smaller, 
oscillations decay rapidly and the echo approaches a con- 
stant that grows with distance. In both cases the depen- 
dence on distance rapidly saturates: the long distance 
regime is reached at d > 4. To the contrary, near the 
critical point (A = 1) the echo decreases logarithmically 
with time (after a short transient), and saturation with 
distance is not attained, which is a signal of long range 
correlations in the environment. 

The strongest decoherence is typically obtained at 
short distances. A simple argument shows that this is 
reasonable for weak coupling: at short times one ex- 
pects an approximately Gaussian decay of the echo of 
the form exp(— a^t^), with a oa g. For d = we have 
~ {2g)^ while in the case of independent environments 
(X 2g^. As a consequence, faster decay of Loo, 11 is ex- 
pected for d = in the weak coupling regime. For the 
echo Loi.io the dependence on distance is opposite to the 
previous case, as shown in Fig. [3) This is expected since 
the Hamiltonians Hqi and Hiq become more different 
as d grows. On the other hand, the saturation with dis- 




FIG. 2: Echo I/oo.ii as a function of time for an Ising chain 
(7 = 1) of A'^ = 100 sites weakly coupled to a two-qubit 
system (g — 0.1). From top to bottom the transverse field 
is A = 0.5, 0.99, and 1.5. The distance between qubits is 
d = (full black), 1 (dashed green), 2 (dotted violet), and 3 
(dash-dotted magenta). In the almost critical case A — 0.99 

we include also d — 4 {■ ■ —, red), and 10 ( ■, orange); 

these curves are not shown in the other plots because they 
are intertwined with the others. For comparison we display 
in full grey the limit of two independent environments. 



tance approaching the limit of independent environments 
is similar to that of Loo,ii(5, t). 

We analyzed the saturation scale in the foUowing way: 
we calculated the echo Loo, 11 for A between 0.1 and 2, 
distances d between 1 and 15, and times tj in the in- 
terval [0, 10]. The squared norm of the vector obtained 
substracting the values of the echo at times tj for dis- 
tance d and for independent environments was taken as 
a measure of the difference between these echoes. A sat- 
uration length I was obtained from the exponential fit 
of the decay of these norms as d approaches the long 
distance limit. The behaviour of Z as a function of the 
external field A is shown in Fig. [D The saturation 
length was found to be related to the correlation length 
^ = |ln(A)|~^ O Our numerical results indicate 

that I ~ 1.1-1-0.21* (0.17 -I- C"i)"2.2^ Indeed, / is clearly 

increased near the critical point, though it does not di- 
verge like ^. The analysis of the echo Loi,io leads to 
analogous results. 

The echoes obtained for the case 7 = 0.1 are qualita- 
tively similar to the ones shown in Figs. [H [3] for the 
Ising chain. There are nevertheless some noticeable dif- 
ferences. These concern mainly the dependence on dis- 



FIG. 3: Echo I/oi,io as a function of time; all parameters are 
the same as in the previous Figure (for d = this echo is 
always equal to 1). 




FIG. 4: Saturation length as a function of the external field 
A (dots), for the Ising chain (7 — 1). The continuous red 
line shows the fit as a function of the correlation length ^ of 
the chain {I ~ 1.1 + 0.21 * (0.17 -I- C"^)"^'^)- These results 
correspond to a chain of 500 spins, and the echoes have been 
smoothed to avoid spurious effects due to oscillations. Be- 
sides, the fact that y > is accounted for by a shift in the 
correlation length, given by ^cff ~ C(-^ + 5-7 10~^). 



tance, which is rather irregular for A < 1. Besides, the 
connection between the saturation distance / and the cor- 
relation length ^ is not evident, even though the satura- 
tion length clearly increases in the proximity of the crit- 
ical point. There are also differences in the shape of the 
decay close to criticality (almost linear for times up to 
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FIG. 5: The echo for two qubits weakly interacting with a 
spin chain of A'^ = 100 sites has revivals or sharp decays due 
to finite-size effects. For a distance d = N/2 we show Loo.ii 
(full black) and I/oi,io (dashed red). When d decreases the 
first peak/decay tends to dissappear (the dotted blue line cor- 
responds to 1/00,11 for d — 30). The result for two qubits cou- 
pled to independent environments is shown for comparison 
(full grey). The parameters are A = 0.99, 7 = 1, <? = 0.1. 



FIG. 6: Envelope of the echo Loo, 11 a function of time for 
two qubits strongly interacting with two different sites of a 
chain with = 100, A = 0.99, 7 = 1, 5 = 50. The plots cor- 
respond to the cases of distance d = (full black), 2 (dotted 
blue), and long-distance limit (full grey). For the case d = 
a part of the fast oscillation is included. Here, in contrast 
to the weak couphng scenario, decoherence is stronger in the 
long distance limit, as the envelope for large d corresponds to 
the square of the one for d = 0. 



t ^ 10), and in the strength of the decoherence process 
(for A < 1 decoherence is stronger than in the case 7=1, 
while for A > 1 it is weaker). 

The previous results refer to times shorter than the 
ones where the finite size of the environment starts to 
play a role. For long times finite size effects become im- 
portant, as shown by the coherence revivals and sharp 
decays of Fig. [5l These imply that the qubits are not in- 
dependent due to environment-induced interactions. By 
decreasing d below N/A the first peak or decay tends to 
disappear, so that the qubits evolve almost independently 
even for long times (provided d is over the saturation 
scale). This will not be the case in the strong coupling 
regime, which will be treated in the next section. 

IV. RESULTS: STRONG COUPLING 

For strong system-environment coupling, the results 
are quite different from the ones above. As observed for 
a single spin system in il5] , the strong coupling regime 
is characterized by an echo with a fast oscillation and a 
slow envelope which for large enough g is independent of 
g. Following the same steps as in [5| we can explain the 
behaviour of the echo Loo, 11 as follows: Consider the evo- 
lution with Hamiltonian Hn = Hqq — g{Zo + Zd). As all 
the frequencies associated with the interaction Hamilto- 
nian are of order g (much larger than typical frequencies 
of the chain), we may approximate: 

e-'"'''\Eo) ~ e'3*(^'>+^'''^e-'*"'\Eo). (20) 

Here H' is the Hamiltonian iJoo reduced to a block di- 
agonal form, with blocks associated with the different 



eigenvalues of the interaction term. The evolution oper- 
ator thus factors in two parts: a fast periodic evolution 
with frequency of order g, and a slow evolution that de- 
termines the envelope of the echo and is governed by an 
effective chain Hamiltonian H' (which does not depend 
on g). This same argument can be adapted to Loi^io. We 
note, however, that this behaviour is a consequence of the 
form of the interaction between system and environment 
and is not generic for strong couplings. For instance, if 
each qubit interacts with a region of the chain with a 
coupling that decreases with distance, the typical result 
is a sudden decay of coherence without oscillations. 

For strong coupling we expect decoherence to increase 
with distance (until saturation) for both echoes. Indeed, 
for d = we have Loo, 11(5, i) = -^o,i(2g, i), with Lo,i 
the single qubit echo. This must be of the same order 
as Lo,i(g,t) because for large enough g the envelope is 
independent of g. On the other hand, in the long dis- 
tance limit Loo, 11 (3, t) — Lq i{g, t), which is smaller than 
La.i{g,t). The comparison between d ^ 1 and d = 0, 
illustrated in Fig. [5] for 7 = 1, A = 0.99, thus leads 
to a result which is quite the opposite of the one ob- 
tained for weak coupling. The distance dependence of 
the echo Loi^io follows a similar pattern, i.e., decoher- 
ence increases with distance, as shown in Fig. [71 

Environment-induced interactions between qubits are 
manifest in Fig. [51 Indeed, for d > appart from the 
fast oscillation with a decaying envelope the echo Loo. 11 
exhibits a beating (or revival) at a time that depends on 
d. This effect is not important for weak coupling and 
can be interpreted as due to the interaction between the 
qubits through the modes of the perturbed chain. On 
the contrary, no revivals are seen for Loi.io- 
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FIG. 7: Envelope of the echo I/oi.io as a function of time 
for two qubits strongly interacting with two different sites 
of a chain with iV = 100, A = 0.99, 7 = 1, 5 = 50. The 
plots correspond to the cases of distance d = 1 (full black), 
2 (dotted blue), and long-distance limit (full grey). For the 
case d = there is no decoherence in this case, while for d > 1 
decoherence turns out to be stronger as distance is increased. 



The difference can be understood by analyzing the ori- 
gin of revivals in terms of the spectrum of the Hamil- 
tonian and the Bogoliubov coefficients, contained in the 
matrices A, g and h in ([T5]). Excitations of the unper- 
turbed Hamiltonian have energies lying between 2|1 — A| 
and 2|1 -I- A|. When a strong external field is applied in 
two sites, two eigenvalues of order g appear; these ex- 
citations are associated to combinations of the original 
fermion operators in the two sites. The remaining ex- 
citations (with eigenvalues of the same order as before) 
can be split in two groups, corresponding roughly to ex- 
citations between and outside the interaction sites (this 
labelling only makes sense because we consider distances 
d N). It turns out that the most populated levels cor- 
respond to the lowest energy excitation (occupying the 
outside region), the excitations in the interaction sites, 
and those lying in the inside region . The high-energy 
excitations are associated to the rapid echo oscillation, 
while the lowest energy excitation has a time scale grow- 
ing with N. The beating in the echo is given by the 
lowest-energy mode in the region between the qubits. 

The situation for Loi.io is different, as seen in Fig. [3 
because for each effective Hamiltonian there is a single 
site perturbation of the chain. In this case the chain 
cannot be split in two regions and the populations of the 
low-energy levels decrease slowly with the energy, in such 
a way that there is not a single frequency associated to 
them. Then, no revivals occur up to times long enough 
for finite-size effects to appear. 

The revivals for the case Loo, 11 appear for A < 1 but 
not for A > 1. This is a consequence of the change in 
the properties of the chain. Indeed, for A < 1 the re- 
vival time is shorter as A is increased because the lowest 
energy associated to excitations between the qubits in- 
creases. For A > 1, the dissappearance of the revivals 



may be related to the fact that many low-energy excita- 
tions become populated, with frequencies comparable to 
the one associated to the beating. 

The first revival time, t^, and the amplitude of the 
peak, Lr, can be analyzed as functions of the distance 
d. We consider first a critical case where Hqq and Hu 
belong to different phases (A = 0.99, 7 = 1, iV = 100). In 
each case, new peaks appear at multiples oitr- For d > 1, 
the revival time grows linearly with distance {tr « 2d), 
while Lr decays as a power law, Lr oc d~^^'^. For dis- 
tances d > N/5 effects associated with the finite size of 
the chain appear altering the regular trend. Examining 
the echo for A far away from the phase transition we 
found a slower decay of the revival peak with d. Fur- 
thermore, in this case tr is not linear with d, but seems 
to grow exponentially. Thus, tr appears to be related 
to transport properties of the chain that are modified by 
varying the external field. This should serve as a warn- 
ing not to picture revivals as the manifestation of a spin 
wave propagating with constant velocity along the chain 
(the energies of the lowest modes between the qubits do 
not generally scale as 1 /d) . 

Changing the anisotropy parameter 7, the envelopes 
display different shapes, heights and revival times. The 
distance dependence of the revival time tr and amplitude 
Lr was studied for A — 0.99 and 7 between 0.1 and 1. We 
found that tr generally grows as a power law. Besides, the 
height of the revival peaks decreases with increasing 7, 
and the dependence of Lr with d maintains a power-law 
decay. We note, however, that the peaks lose definition 
for large d and small 7. 



V. CONCLUSIONS 

We studied the evolution of a system of two qubits in- 
teracting locally with an XY chain. Decoherence and dis- 
entanglement are determined by a Loschmidt echo, which 
we computed using previous results 0, [l3| • The formulas 
we obtained for the echoes enabled us to study a family 
of initial states including the four Bell states. We focused 
our analysis on the dependence of the echo with the dis- 
tance between the interaction sites. The cases we studied 
show a rich variety of results. For instance, for strong 
coupling decoherence increases with increasing distance 
for all Bell states. On the contrary, for weak coupling 
decoherence may increase or decrease with distance ac- 
cording to the initial state chosen. 

In the regime of weak system-environment coupling 
the dependence of the echo with distance saturates very 
fast except for the critical case. The saturation length 
I characterizing the approach to the long distance limit 
seems to depend on the correlation length of the chain ^. 
This result is interesting as it relates an equilibrium quan- 
tity (the correlation length ^) with a dynamical quantity 
I. However, this relation was only found for the Ising 
case, and will be further analyzed elsewhere. 

The strong coupling regime is characterized by fast os- 
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dilations with a slowly decaying envelope. For some ini- 
tial states, this envelope shows a beating which can be 
interpreted as an interaction of the qubits through the 
chain. Its time scale is determined by the modes of the 
chain in the region between the qubits, which depend on 
the distance between them and on the values of the pa- 
rameters in the chain Hamiltonian (indeed, the beatings 
were found only for A < 1). We note that in this regime 
the fast oscillation of the echo continually provokes de- 
cays and revivals of the entanglement between qubits. 
This is just a dynamical transfer of this entanglement 
back and forth from the system to its immediate vicin- 
ity with a frequency given by the coupling constant g. 
The true loss of entanglement is produced by the decay 
of the oscillation. If the initial state is mixed, this decay 
can lead to true sudden death of entanglement. When 
the distance between qubits is short, the sudden death 
may be followed by a sudden revival due to the beat- 
ing in the echo. For very long distances, entanglement 
loss becomes irreversible, unless we consider times which 
are long enough for the finite size of the environment to 
become manifest. 

Finally, it is worth pointing out that there are formal 
analogies between the model we solved and other deco- 
herence models. In fact, decoherence for two qubits at 
distance d in an initial state of the form ajOO) + (3\11) 
is identical to that of a single qubit interacting non- 
locally with sites and d (which was partially studied 
in [^). On the other hand, for the case of states of 
the form a|01) + /3|10), the relevant echo Loi,io can be 
mapped onto the one of another equivalent problem: a 
spinless particle occupying discrete positions along the 
chain (with position eigenstates \xj), j = 0, . . . , iV — 1). 
If this particle modifies the effective magnetic field for 
the spin in the site it occupies, cat states of the form 
{\xj) + |.Tj+(i))/-\/2 will decohere as a consequence of the 
interaction. The resulting reduction of the off-diagonal 
terms will be precisely given by expression l|17p . 

Note added: After the submission of this paper, we 
became aware of the work in [23|, which shows how to 
calculate efficiently the remaining elements of the density 
matrix (/9oi,ii)- 



VI. APPENDIX A 

In this section we show how the formula (fTS]) for the 
echo can be obtained. Taking into account the relation 
between the operators that diagonalize the different 
effective Hamiltonians, 



(0) 



E 

k 



(1) 

ijkVk 



the two vacuum states \Eq) (for the unperturbed Hamil- 
tonian) and \E'r) (for the perturbed one) can be con- 
nected by [ll|: 

|i?o)cxcxp|i7f«G^^w||i?^) (21) 



with G = — g~"'^h (for the sake of simplicity, we shall first 
assume that g is invertible, and sketch the most general 
case afterwards). The echo can then be calculated from: 



L{t) = \{E^\e-^"-'\E^)\' 

cx \{E'Q\e-^-^^'^e''-''^'^^e^^'^^' \E'^) p 



(22) 



where in the last expression A is a diagonal matrix con- 
taining the energies corresponding to the different par- 
ticles, and all super-indices are omitted as all operators 
and matrices refer to the perturbed Hamiltonian Hi . By 
introducing two identities in terms of fermionic coher- 
ent states between the exponentials and integrating two 
times, [l9| we obtain: 



L{t) cx 



dai . . . dajsfdPi . . . d[3 



N 



exp <; iaG*a -t- fie'^^a - ^^GP 



(23) 



where c5, (3 are Grassman A^-tuples. This is a Gaussian 
integral, that can be solved to: 



det I ^j4A _Q 



(24) 



Using properties of the determinant and the fact that 
A is diagonal, and imposing that L(0) = 1 the echo can 
be rewritten as: 



ioi,io(i) = |det(ge^*V+he-^*^ht)| 



(25) 



The formula can be simplified further: since g, h are real 
we have that g±h is orthogonal, and using this we obtain 
the final expression ^ 



L{t) = |det(g-Hhe*^*)| 

In case g is not invertible, the relation pT|) between 
the vacuum states can be g eneralized by using interme- 
diate sets of operators [l3|. By the singular value de- 
composition, g = UDV with D diagonal, U, V orthog- 



onal. We define new fermionic operators ^ 



— Vff ^h The linear transformation between the 
^t(«) j^QTf^ i^as £) instead of g, and the vacuum states are 
the same because the transformation does not mix cre- 
ation and annihilation operators. We can assume Dj = 
for j < ji, and the remaining eigenvalues to be nonzero. 



By interchanging particles with holes (^^'^■' ^ Cj^"') for 
every index such that Dj = we obtain a linear trans- 
formation with an invertible matrix. The calculation of 
the echo follows the same steps as before, except that it 
is necessary to treat indices j < ji separately. After the 
Gaussian integration, we are left with an expression of 
the form (|24p in which only indices j > ji appear. The 
desired result can be achieved by conveniently introduc- 
ing some rows and columns in the matrix, in such a way 



^t(0)^ 



8 



to include the parts of the matrices with j < ji without 
changing the value of the determinant. 

The present derivation cannot be easily extended to 
more general cases, with two different evolution operators 
as in ([7]) , for in that case there are three sets of operators 
that must be related. The steps taken in the calculation 
above then lead to a result that involves a sum over an 
exponential number of terms, and so cannot be efficiently 
evaluated. 



VII. APPENDIX B 

Here we sketch the derivation of the formula (|17p for 
the decay of the off-diagonal terms poi.io- We start from 
expression (llOp . which involves two sets of particle op- 
erators: one which we note as ri^^\ diagonalizing the 
effective Hamiltonian Hiq, and the set 77'^"^ for the un- 
perturbed chain Hamiltonian Hqq. These two sets are 
connected by a real Bogoliubov transformation. 

The translation operator can then be written as [2l| : 



(26) 



so that the echo is given by: 

LoiMt) = (^ole"^ \Eo) ' (27) 



with 



g^-ffiotj^Wg-i-f/iot^ 



(28) 



This leads us to an echo of exactly the same form as that 
in Appendix A, except that the energy matrix A is re- 
placed by a matrix with the eigenvalues of the translation 
operator, and the set of particles 77^^^ is replaced by the 
Q!fc. The main difference here is given by the fact that 
the Bogoliubov transformation between the ak and the 
77^*^^ has time-dependent complex coefficients: 



(0) 



k 



h'jkal 



with 



ge**^g* + he 



h' 



ge**^h* + he"^*'^-' 



-**^h*, 



(29) 



(30) 
(31) 



Following steps similar to those in Appendix A, we finally 
obtain for the echo the expression ([TT]): 



(t) = det(g'r-V'+h'T'^h") 



^01,101 



where the diagonal matrix t contains the eigenvalues of 
the translation operator T; the structure is the same as 
in In this way the decay of the off-diagonal terms 

can once more be written as the determinant of an TV x TV 
matrix. This expression is not as simple as the final one 
obtained in Appendix A; the reason for this is that the 
matrices g±h were orthogonal, while the matrices g'±h' 
do not satisfy a similar condition. 
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